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A theory of three-dimensional (3D) hypothetical magnonic crystal (conceived as the magnetic 
counterpart of the well-known photonic crystal) is developed and applied to explain the existence of 
a spin-wave frequency gap recently revealed in low-doped manganites Lai- x Ca x Mn03 by neutron 
scattering. A successful confrontation with the experimental results allows us to formulate a working 
hypothesis that certain manganites could be regarded as 3D magnonic crystals existing in nature. 

PACS numbers: 75.30.Ds, 75.25.+Z, 
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INTRODUCTION 



Though the first study of electromagnetic wave prop- 
agation in periodic structures, written by lord Rayleigh, 
was published already in 1887, not until very recently 
have photonic composites suddenly raised an extremely 
keen interest. The research in this field was triggered 
by the studies by Yablonovitch and John Q, 0, pub- 
lished in 1987 and anticipating the existence of complete 
energy gaps in electromagnetic wave spectra in three- 
dimensional periodic composites, henceforth referred to 
as photonic crystals. These can be used for fabricating 
new optoelectronic devices with photons acting as trans- 
port medium. The so-called left-handed materials, an 
example of periodic structures characterized by negative 
effective refractivity Q , demonstrate how much the prop- 
erties of this kind of structure can differ from those of 
homogeneous materials. Another type of periodic com- 
posites are structures composed of materials with dif- 
ferent elastic properties; showing an energy gap in their 
elastic wave spectrum, such composites are referred to 
as phononic crystals |4|. Recently, attention has been 
focused on the search of photonic and phononic crystals 
in which both the position and the width of the energy 
gap could be controlled by external factors, such as ap- 
plied voltage or magnetic field. Attempts are made to 
create photonic crystals in which one of the component 
materials would be a magnetic [3, |a| . 

A magnetic periodic composite, conceived as the mag- 
netic counterpart of a photonic crystal, consists of at 
least two magnetic materials, with magnon acting as in- 
formation carrier, therefore such periodic magnetic com- 
posites can be referred to as magnonic crystals (MC). 
Studies of 2D magnonic crystals have already been re- 
ported H 0, 0, 03, E3] with scattering centers in the 
form of " infinitely" long cylinders disposed in square lat- 
tice nodes (cylinder and matrix materials being two dif- 
ferent ferromagnetics) and spin wave spectra showing the 
anticipated gaps. Further research was focused on mag- 
netic multilayer systems, which can be regarded as ID 
magnonic crystals |l3l Il ll| . In this paper, we present nu- 



merically calculated - and, to our best knowledge, not 
yet reported in the literature - band structures of three- 
dimensional magnonic crystals. Due to the complexity 
of the problem, only the simplest model of 3D magnonic 
crystal is considered here, represented by a system of 
ferromagnetic spheres (which act as scattering centers) 
disposed in the nodes of a cubic bcc crystal lattice and 
embedded in a different ferromagnetic material (matrix). 
Both the exchange and dipolar interactions are taken into 
account in our calculations, which are based on the plane 
wave method and use the linear approximation. As a 
conclusion, we propose a new magnonic interpretation of 
experimental results obtained through neutron scattering 
on spin waves in doped manganites [TrL Hfil Il7j . 



THEORY OF 3D MAGNONIC BAND 
STRUCTURE 



a) 




FIG. 1: a) The 3D periodic structure studied in this paper; 
the structure consists of ferromagnetic material A spheres 
embedded in a material B matrix (materials A and B having 
different magnetic properties) and disposed in the nodes of 
a bcc lattice, (b) The Brillouin zone corresponding to the 
considered structure, and the path (PFH N P, highlighted) 
along which the magnonic band spectra are calculated. 



Let's consider an ideal periodic structure consisting of 
spheres of ferromagnetic A embedded in a matrix of fer- 
romagnetic B. The spheres are assumed to form a 3D 
periodical lattice of bcc (Fig. la) type. A static mag- 
netic field, Ho, is applied to the composite along the 
z axis, and assumed to be strong enough to saturate 
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the magnetization of both materials. The lattice con- 
stant is denoted by a; the filling fraction, / = | ir R 3 a~ , 
is defined as the volume proportion of material A in a 
unit cell. Ferromagnetics A and B are characterized by 
two material parameters: the spontaneous magnetization 
(Ms,a and Ms,b), and the exchange constant (Aa and 
Ab); both these parameters depend on the position vec- 
tor r = [x, y, z): 



M s {r) = M StB + {M s , A -M StB )S{f), 
A(r) = A B + (A A -A B )S(r), 



(1) 



the value of function S(r) being 1 inside a sphere, and 
beyond. 

In the classical approximation, spin waves are de- 
scribed by the Landau-Lifshitz (LL) equation, taking the 
following form in the case of magnetic composites (with 
damping neglected): 



dM{f,t) 

m 



= ytx M(r,t) x H ef f(r,t), 



(2) 



where magnetization M(r, t) is a function of position vec- 
tor r and time t; H e f -t (r , t) stands for the effective mag- 
netic field pllauillii acting on magnetization M(r,t): 



H eff (f,t)=H z + h(f,t) + 
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)M(r,t):(3) 



z is the unit vector along the z axis; h(f, t) is the dy- 
namic magnetic field resulting from dipolar interactions; 
the third component represents the exchange field. The 
magnetization vector can be represented as the sum of 
its static and dynamic components; the former, M$z, 



J 



is parallel to the applied field; the latter, m(r, t) — 
rh(r) exp(— iu>t), ui denoting the precession circular fre- 
quency, lies in the plane (x,y): 



M(f,t) = M s z + m{r,t). 



(4) 



The dynamic dipolar field, h, must satisfy the magneto- 
static Maxwell equations: 



V x h{r) = 0, 
V • (h(r) + m(f)) = 0. 



(5) 



In magnonic crystals, the position-dependent coefficients 
in i.e. Ms and A, are periodic functions of the po- 
sition vector, which allows us to use in the procedure 
of solving the LL equation the plane wave method, 
described in detail in our earlier papers dGii (dealing 
with two-dimensional magnonic crystals). Following this 
scheme, we proceed to Fourier-expanding all the periodic 
functions of the position vector, i.e. the spontaneous 
magnetization, Ms, and parameter Q defined as follows: 



Q 



2A 



(6) 



The dynamic components of the magnetization can be ex- 
pressed as the product of the periodic envelope function 
and the Bloch factor, exp(ikr) (k denoting a 3D wave 
vector); the envelope function can be transformed into 
the reciprocal space as well. Including all the expansions 
into (0) and J5j leads to the following infinite system of 
linear equations for Fourier coefficients of the dynamic 
magnetization components, m x ^(G) and m y ^(G): 



i Slm x% {G) = m y% {G) + Y J 



(k y + G' y )(k x + G' x )m x% {G>) + [k y + G' y fm y ~ k {G') 



M S {G-G') 



Ho\k + G'\ 2 

XEP + G') ■ (k + G") — (G — G") ■ (G - G')]M S {G - G")Q{G" - G')m y% (G'), 



G' G" 



_ _ (k v + G')(k x + G')m AG') + (k x + G') 2 m AG') ^ ^ 

iflm AG) = -m AG) - V — v — x> yk S L — — — -Ms(G - G ) 

yk\ > X kV > Z_, H \k + G'\ 2 K ' 

G 

~J2J2^+ • (* + - - • @ - G')Ws(G - G")Q(G" - G')m xi {G')- 

G' G" 



(7) 



k x , k y and G x , G y denoting the Cartesian components of 
the wave vector k and a reciprocal lattice vector G, re- 
spectively; a new quantity introduced in (JJJ is f2, hence- 



forth referred to as reduced frequency: 



h\lM)Ho' 



(8) 
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The Fourier coefficients of spontaneous magnetization 
Ms and parameter Q are calculated from the inverse 
Fourier transformation; in the case of spheres, the re- 
sulting formulae for Mg read as follows: 

f M S}A f + M StB (l - f), forG = 

M (C\ ) f(M M MMGR)-{GR)cqz{GR)] 
M S {G) = < f{M s ,A - Ms,b) [GR) 3 

[ for G>0 

and similarly for Q; R is the sphere radius (Fig. QJi). 

Obviously, the numerical calculations performed on the 
basis of J7J involve a finite number of reciprocal lattice 
vectors G in the Fourier expansions; however, we have 
made sure that the number used is large enough to guar- 
antee good convergence of the numerical results. As in- 
dicated by an analysis performed, a satisfactory conver- 
gence is obtained already with 343 reciprocal lattice vec- 
tors used. 

NUMERICAL RESULTS AND CONFRONTATION 
WITH THE EXPERIMENT IN MANGANITES 



a) b) 
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FIG. 2: Magnonic bands computed on the basis of J7[ for 
a hypothetic periodic arrangement of ferromagnetic droplets 
embedded in different ferromagnetic medium; the sphere- 
shaped droplets are assumed to form a regular crystal lattice 
of bcc type, (a) Magnonic branches plotted versus the lattice 
constant, a, of the droplet bcc structure; the filling fraction 
is fixed at / = 0.33. A frequency gap (the shaded region) is 
found to exist, and move down with increasing a. The super- 
imposed empty rectangle represents the gap revealed experi- 
mentally in doped manganites |Ts||. (b) The magnonic energy 
branches (in the same crystal) plotted versus the filling frac- 
tion, with the lattice constant fixed at a — 42A. 



An energy gap in spin wave spectra has re- 
cently been found experimentally in doped manganites 
La\- x C 'a x MnOz and La\- x Sr x MnO^ at low (Ca or Sr) 
ion concentrations (2% < x < 10%) [HE1]. The ap- 
pearance of this gap was suggested to be due to doping, 
which results in a two-phase magnetic structure, with fer- 
romagnetic (F) droplets (associated with the Mn A+ ions) 
embedded in a canted antiferromagnetic (CAF) medium 



(showing a non-zero magnetic moment along the direc- 
tion of the droplet magnetic moments) . Hence, the two 
separate energy branches found in the neutron scatter- 
ing experiments are interpreted as corresponding to two 
independent spin-wave excitations, each propagating in 
one of the two phases. In this paper, we propose an 
alternative magnonic explanation of the spin-wave spec- 
trum energy gap found in these experiments. We antici- 
pate that the dispersion branches obtained in the exper- 
iment constitute the bottom part of the collective spin- 
wave spectrum, involving spin waves propagating con- 
currently through both phases of a magnetic composite, 
resulting from a regular arrangement of ferromagnetic 
droplets (material A) in a matrix with different mag- 
netic properties (material B). Then, the existence of a 
spin-wave spectrum gap becomes a natural consequence 
of the structure periodicity, as in the case of photonic 
or phononic crystals. Below we present results of our 
calculations (based on the method discussed in the pre- 
ceding paragraph) of the spin- wave spectrum in a droplet 
magnonic crystal, the necessary material parameter val- 
ues being as in the experimental studies 
results are compared to the spectrum obtained experi- 
mentally in Lao.9Cao.1A/nO3. 

Let's assume that the spherical ferromagnetic droplets 
form a bcc lattice (Fig. HJi), which we believe to be a 
good approximation of the actual droplet structure. The 
spin-wave branches will be computed along the follow- 
ing path in the irreducible part of the first Brillouin 
zone: P = 7r/o(l,l,l) -> V = vr/a(0,0,0) -> H = 
w/a[2, 0, 0) -> N = 7r/o(l, 1, 0) -» P = n/a[l, 1, 1) (see 
Fig. m>). The necessary material parameter values, i.e. 
those of the spontaneous magnetization and the exchange 
constant, in both magnetic phases of our model magnonic 
crystal are as in the experimental studies referred to 
above. The paper by Hennion et al. 0] provides in- 
formation on the difference, Am, between the values of 
magnetization along vector c in the droplets and in the 
matrix: Am = 0.7 ± 0.2/ib- Assuming that the mag- 
netic moment in droplets comes mainly from Mn A+ ions 
with spin S — 3/2, we get droplet spontaneous magne- 
tization value Ms, a = 0.46 10 6 Am -1 ; consequently, the 
matrix magnetization [Ms, a l ess 0.7/x# = 0.22Am _1 ) is 
M S B = 0.24 ItfAmT 1 . The following lattice constant 
values in Lao.gCao iMnOs are assumed: a = 5.465A, 
b = 5. 621 A , and c = 7.725A [13 ■ The exchange con- 
stant values in droplets and in the matrix, A a and Ab, 
respectively, are determined on the basis of the exchange 
integral, Ji(F) = ImeV, and the stiffness constant, 
D w IbmeVA 2 , as estimated in the study by Biotteau 
et al. 16]; the resulting droplet exchange constant value 
is A A = S 2 J 1 {F)/a c [meVA- 1 ] = 0.09 • 10~ n JmT 1 , a c 
being the lattice constant recalculated for the manganite 
pseudocubic elementary cell. The matrix exchange con- 
stant value, Ab, is determined by fitting the weighted 
mean value of the droplet and matrix stiffness constants 
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to the stiffness constant value estimated in Jf6|. We 
assume As = 0.012 JmT x , which, at filling fraction 
/ = 0.2, corresponds to mean value D = 18.31 meVA 2 , 
close to the value estimated by fitting to the experimen- 
tal results at temperature f5°K. The assumed average 
value of the isotropic internal field, fioH = OAT, lies be- 
tween the values found by Savosta et al. H3 (0.064T) 
and Yates et al. (0.35T). 

a) b) 
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FIG. 3: (a) The (a, /) plane with two sets of lattice con- 
stant a and filling fraction / values in a magnonic crystal: 
one, represented by the solid line, corresponds to an energy 
gap whose width is exactly Aui — 0.2meV; the other, rep- 
resented by the dotted line, corresponds to a gap centered 
around loq = 1.65meV. The cross-section points indicate the 
a and / values at which the computed gap fits the experimen- 
tal one in terms of both width (Au) and position (u>o)- (b) 
Magnonic bands calculated on the basis of our model across 
the Brillouin zone, the assumed lattice constant and filling 
fraction values being a = 38. 5 A, and / = 0.13, respectively. 
The gap which is found to exist has position and width exactly 
fitting the experimental data reported in 

Figure |2t shows the computed magnonic band struc- 
ture plotted against the droplet lattice constant. The 
gap center (u>o) is found to descend, and the gap itself 
(the dark gray region) to narrow down as the lattice con- 
stant increases. The superimposed rectangle represents 
the 'experimental' gap reported in ^(|). The gap cen- 
ter found experimentally corresponds to our results at 
a = 42A; however, the corresponding gap width result- 
ing from our computations is much above the experi- 
mental result, which implies the necessity of fitting the 
filling fraction as well. Fig. |2h shows the magnonic band 
structure plotted against the filling fraction / (the lat- 
tice constant being fixed at a = 42 A). The gap is found 
to open at / = 0.10 to reach its maximum width at 
/ « 0.42 and vanish at f — 0.64, which means that one 
gap width value corresponds to two different filling frac- 
tion values. Note also that the gap center moves up as 
the filling fraction increases. Hence, the a and / val- 
ues corresponding to an energy gap that fits the experi- 
mental results will be determined in two steps (see Fig. 
|3Ji): first, we shall find the set of a and / values corre- 
sponding to the experimental gap width, Aui = 0.2meV 
(the solid line in Fig|3^); next, another line (the dot- 
ted one in the Figure) is plotted, corresponding to the 



gap centered at the experimental value, ujq = 1.65meV. 
The cross section of these two lines indicates the a and 
/ values at which both ujq and Aoj are fitted. As the 
lines cross at two points: (a = 38. 5A, f = 0.13) and 
(a = 56. 5A, f = 0.61), two solutions are obtained, 
one corresponding to the droplet radius R = 9.6A, the 
other to R = 23. 6A. The former solution proves close 
beyond expectation to the experimental result reported 
in 0| (Fig. 9), where La . 9 Ca .iMnO3 droplet radius 
is estimated at R — 9.37A The corresponding mean 
droplet spacing estimated in is d m f» 36A, also quite 
close to our result, a = 38. 5A. Therefore, the complete 
set of structural and material parameters characteriz- 
ing our droplet magnonic model of Lao.gCao,\Mn03 is: 
M S , A = 0.46 lOMm" 1 , M s .b = 0.24 10 6 Ato-\ A a = 
0.09 10" 11 Jm" 1 , A B = 0.012 10" 11 Jm" 1 , hqH = OAT, 
a = 38. 5A, f = 0.13. The corresponding full spin-wave 
spectrum is shown in Fig. |3t>- An energy gap, existing 
indeed between the first band and the second one, cor- 
responds very well, in terms of both position and width, 
to that found experimentally 0] . 

This result allows to propose a working hypothesis that 
low-doped manganites could be regarded as magnonic 
crystals existing in nature. Obviously, further research is 
required for verification of this hypothesis; therefore, in 
another paper, we are going to investigate the effect of 
non-spherical droplet shape on the magnonic spectrum, 
and to consider a real canted antiferromagnetic matrix, 
rather than the effective ferromagnetic one considered 
here. However, already the results obtained here on the 
basis of the approximate model are extremely promising. 



Acknowledgements 

The present work was supported by the Polish State 
Committee for Scientific Research through projects KBN 
- 2P03B 120 23 and PBZ-KBN-044/P03-2001. 



[1] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987). 
[2] S. John, Phys. Rev. Lett. 58, 2486 (1987). 
[3] J. B. Pendry, D. R. Smith, Physics Today 57, 37 (2004). 
[4] M. S. Kushwaha, P. Halevi, L. Dobrzynski, B. Djafari- 

Rouhani, Phys. Rev. Lett. 71, 2022 (1993). 
[5] A. Saib, D. Vanhoenacker- Janvier, I. Huynen, A. Enci- 

nas, L. Piraux, E. Ferain, R. Legras, Appl. Phys. Lett. 

83, 2378 (2003). 
[6] S.A. Nikitov, Ph. Tailhades, Optics Communications 

199, 389 (2001). 
[7] L.D. Landau and E.M. Lifshitz, Phys. Z. Sowjetunion 8, 

153 (1935). 

[8] J. O. Vasseur, L. Dobrzynski, B. Djafari-Rouhani, H. 

Puszkarski, Phys. Rev. B 54, 1043 (1996). 
[9] M. Krawczyk, H. Puszkarski, Acta Physica Polonica A 

93, 805 (1998). 



■5 



[10] H. Puszkarski, M. Krawczyk, Solid State Phenomena 94, 
125 (2003). 

[11] Y.V. Gulyaev, S.A. Nikitov, L.V. Zhivotovskii, A. A 
Klimov, P. Tailhades, L. Presmanes, C. Bonningue, C.S. 
Tsai, S.L. Vysotskii, Y.A.Filimonov, J.E.T.P. Letters 77, 
670 (2003). 

[12] M. Shamonin, A. Snarskii and M. Zhenirovskiyy, 

NDT&E International 37, 35 (2004). 
[13] H. Puszkarski, M. Krawczyk, Physics Letters A 282, 106 

(2001). 

[14] V.V. Kruglyak, A.N. Kuchko, J. Magn. Magn. Mater. 

272-276, 302 (2004). 
[15] M. Hennion, F. Moussa, G. Biotteau, J. Rodriguez- 



Carvajal, L. Pinsard, A. Revcolevschi, Phys. Rev. Lett. 

81, 1957 (1998). 
[16] G. Biotteau, M. Hennion, F. Moussa, J. Rodriguez- 

Carvajal, L. Pinsard, A. Revcolevschi, Y. M. Mukovskii, 

D. Shulyatev, Phys. Rev. B 64, 104421 (2001). 
[17] P. Kober-Lehouelleur, F. Moussa, M. Hennion, A. 

Ivanov, L. Pinsard-Gaudart and A. Revcolevschi, Phys. 

Rev. B 70, 144409 (2004). 
[18] M.M. Savosta, P. Novak, Phys. Rev. Lett. 87, 137204 

(2001). 

[19] K.A. Yates, C. Kapusta, P.C. Riedi, L. Ghivelder, L.F. 
Cohen, J. Magn. Magn. Mater. 260, 105 (2003). 



